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ABSTRACT 



Context. Whether or not there exists a physical upper mass limit for star clusters is as yet unclear. For small cluster samples the mass function 
may not be sampled all the way to the truncation, if there is one. Data for the rich cluster population in the interacting galaxy M51 enables us 
to investigate this in more detail. 

Aims. Using HST/ACS data, we investigate whether the cluster luminosity function (LF) in M5 1 shows evidence for an upper limit to the mass 
function. The variations of the cluster luminosity function parameters with position on the disk are addressed. 

Methods. We determine the cluster LF for all clusters in M5 1 falling within our selection criteria, as well as for several subsets of the sample. 
In that way we can determine the properties of the cluster population as a function of galactocentric distance and background intensity. By 
comparing observed and simulated LFs we can constrain the underlying cluster initial mass function and/or cluster disruption parameters. A 
physical upper mass limit for star clusters will appear as a bend dividing two power law parts in the LF, if the cluster sample is large enough to 
sample the full range of cluster masses. The location of the bend in the LF is indicative of the value of the upper mass limit. The slopes of the 
power laws are an interplay between upper mass limits, disruption times and evolutionary fading. 

Results. The LF of the cluster population of M5 1 is better described by a double power law than by a single power law. We show that 
the cluster initial mass function is likely to be truncated at the high mass end. We conclude from the variation of the LF parameters with 
galactocentric distance that both the upper mass limit and the cluster disruption parameters are likely to be a function of position in the galactic 
disk. At higher galactocentric distances the maximum mass is lower, cluster disruption slower, or both. 



Key words, galaxies: star clusters - galaxies: individual: M51 - star clusters: general 



1. Introduction 

The distinction between star clusters and galaxies used to be 
fairly clear, but recently the gap between both has been slowly 
filled up with objects in the intermediate mass range. The most 
massive star clusters were the Milky Way's old globular clus- 
ters, with masses going up to a few 10 6 M G . The least massive 
galaxies known were galaxies of the type of Fornax, having 
masses of about 10 7 M . Many objects are now determined to 
be in between these distinct mass ranges, e.g. 'Ultra Compact 
Dwarfs' (e.g. iMieske et all 12004 Iffiiker et al.lll999l) . 'stellar 
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superclusters' (e.g. Hil ker et al. 1999) and 'Dwarf-Globular 
Transition Objects' (e.g. lHasegan et al.ll2.Q05b . Whether these 
are all different names for the same objects or not remains to 
be seen. 

Are these new classes of objects really formed in distinct 
formation processes, or is there rather a continuous transi- 
tion between all, a nd are star clusters just a more common 
mode of formation? lKissler-Patig et aL ( 20061) put young mas- 
sive star clusters (M — 1O 65 ~ 7 M ) on the scaling relations 
(relations between mass, velocity dispersion and radius) of 
other hot stellar systems and find that these clusters do fol- 
low scaling relations derived for the least massive galaxies, i.e. 
dwarf elliptical galaxies. Less massive cl usters do not show a 
mass-radius re lation. Si mulations by e.g. Fellh auer & Kro upa 
d2005l) and iBekki et all d2004 support the idea that the most 
massive cluster-like objects may have formed from merging 
smaller clusters, thereby establishing a mass-radius relation 
dKissler-Patig et alj|2006l ). 
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As this would imply that the formation scenarios of sys- 
tems considerably more massive than ~ 10 6 M are really dis- 
tinct from 'normal' star clusters, it is interesting to see whether 
we can also find evidence for a maximu m mass for the star 
clust er-like objects. As pointed out by e.g. ljordan et al. ([2007) 
and lWaters et alJ ([2006, and references therein), the universal- 
ity of the globular cluster luminosity function (GCLF) of glob- 
ular clusters in other galaxies is much easier to achieve when 
including a fundamental upper mass limit, in the form of an 
exponential cut-off, for star clusters in dynamical cluster evo- 
lution models. This also holds for the globular cluster system 
of the Milky Way fe.g. lBurkert & SmithlExjoT) . The universal 
turn-over in the globular cluster mass function is achieved more 
easily with cluster evolution models when including a exponen- 
tial cut off in the mass function around 10 6 M . 

To obtain detailed information (for example the mass, age 
and extinction) of a single cluster of which the stars are not re- 
solved requires photometry in a wide range of broadband filters 
(UV to IR, e.g. Ide Grijs et alJ (12005b . and references therein), 
or spectroscopy. In the study of star cluster populations, where 
spectroscopy is a very time consuming process, the cluster 
LF is a useful tool. This is especially true when photometry 
is only available in a small number of passbands and obtain- 
ing the mass, metallicity and extinction of separate clusters is 
hampered by the age-metallicity-extinction degeneracy. In such 
cases, the LF can be used as a statistical tool, giving only infor- 
mation about the population as a whole. 

In various environments the LF can be approximated well 
by a power law distribution (N dL oc L~"dL ), with a slope 
(-a) between -1.8 and -2.4 (e.g. lLarsenl (120021) : Ide Griis et al. 
l2003h . The translation from the LF to the cluster initial mass 



function (cluster IMF) is non-trivial, however, especially when 
clusters spanning a range of ages are present. D uring its life 



time a cluster fades due to stellar evolution (e.g. ISchulz et al 
20021) an d loses stars due to a vari e ty of disruption ef - 



fects (e.g. Baumgardt & Makinol (12003b : lLamers et alJ d2005l) : 
Gieles et al] d2006ci l2007aF T Nevertheless using assumptions 
for the cluster formation rate, the cluster IMF and cluster dis- 
ruption parameters, one can construct expe cted LFs and subse - 
quently compare these with observed LFs (IGieles et al.ll2005l) . 

In three galaxies the LF of the rich star cluster systems ap- 
pears to be better described by a double power law, i.e. two dis- 
tinct parts, both described by a power law, which are separated 
by a bend at some absolute magnitude. This bend magnitude 
differs from galaxy to galaxy. Whitm ore et al. (1999) found for 
the "Antennae" (NGC 4038/4039) a bend at M v =* -10 (stan- 
dard UBVI), with on the faint side a shallo wer slope (~ -1.7) 
than on the bright side (~ -2.7). For M51. lGieles et alJ d2006ab 
found hints for a double power law behavior of the LF for a sub- 
sample of clus ters in the inner part of the disk, which were later 
confirmed by lGieles et a l. (2006b) for a 5 times larger sample, 
based on the same data as used in this s tudy. The slopes on bot h 



sides are similar to the slopes found by Whit more et al. ( 1999) 



but the bend occurs about 1.6 magnitudes fainter. With the dis- 



tance m odulus of the Antennae newly derived by lSaviane et al 
(2008), whi ch differs by 1 .7 magnitudes fro m the older deter- 



0. 1 magnitude. In Giele s et al. (2006a) it is shown that the LF 
of the star cluster population of NGC 6946 is also better fit with 
a double power law, with parameters comparable to M5 1 . Note 
that the slopes at the faint end of the LF for all these galaxies are 
similar to the slopes found for populations with a single power 
law distribution. If the observed or formed number of clusters 
is not high enough to sample the LF up to the point where the 
bend occurs, the bend will not be detected. The galaxies for 
which a double power law LF for their clusters has been re- 
ported are all typified by a high star formation rate. 

Whereas the bend in the LF of the "Antennae" was origi- 
nally Jntergretedas being caused by a bend in the mass func- 



tion, Zha ng & Falll (119991) already mentioned a truncated mass 
function as a possible explan ation for the double power law 
behavior. IGieles et al. (2006a) have shown with analytic clus- 
ter population models that such a bend can occur if the high- 
est occurring cluster mass is no longer determined by the size 
of the sample (i.e. the highest mass present in a cluster popu- 
lation is determined statistically). For e xample, in the case of 



the LMC and SMC (Hunter et al. 2003), the cluster IMF does 



not show signs of an upper mass limit, and therefore sampling 
this mass function results in a maximum mass in the sample 
that is determined by the number of clusters and the slope of 
th e power law cluster IM F. This is argued to be generally true 
bv lWeidner et al.1 (|2004), who argue that the maximum cluster 
mass is a function of the star formation rate only. On the other 
hand, if there exists a physical upper mass limit for star clus- 
ters and the cluster IMF is sampled up to this limit the LF will 
appear in a sh ape more closely resembling a double power law 
dLarsenll2006l) . 

If the cluster IMF is truncated, the physical reason for the 
maximum mass may have its imprints on the LF of subsets of 
the population. Variations of the maximum mass with galac- 
tocentric distance or surrounding ambient densities can hold 
important clues about the formation and evolution of star clus- 
ters. 

Using deep observations by HST using the Advanced 
Camera for Surveys (ACS) of the interacting, face-on, spiral 
galaxy M51, which cover a large fraction of the disk of the 
galaxy and have a resolution of / /05, we are able to selec t 
a sample of thousands of clusters (iScheepmaker et al.l 120071) . 
Those clusters range in luminosity from the magnitude of sin- 
gle bright stars (distinguished on the basis of their spatial ex- 
tent) to the brightest clusters found in spiral galaxies. The size 
of this sample gives us the opportunity to divide the clusters in 
smaller subsets and still have enough statistics for the obtained 
LFs. In this way we are able to study trends with e.g. galacto- 
centric distance and the spiral arm structure of the ga laxy. The 
same observations are used by iHwang & Led d2008l) to deter- 
mine the LF of luminous clusters. They find a single power law, 
but only i nclude clusters brigh ter t han the bend magnitu de, as 
found by IGieles et al. (2006a) and IGieles et al. I d2006bh . with 
a power law slope which agrees with the bright end slope of 
Gieles et alJ d2006allbh . 



mination by Whit more & Schweizerl d 19951) . the difference in 
bend magnitudes between M51 and the Antennae reduces to 



In § |2] we describe the data, source selection and photom- 
etry. The luminosity function of clusters in M51 will be pre- 
sented in § [3j for the whole population as well as for certain 
subsets. In § [4] we perform Monte Carlo simulations of star 
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Table 1. The data set of the observations of M51, taken under 
proposal ID 10452. 



Filter 


Exposure time (s) 


Limiting magnitude 


F435W (B) 


4 x 680 = 2720 


27.3 mg 


F555W (V) 


4x340 = 1360 


26.5 my 


F814W(I) 


4x340 = 1360 


25.8 m T 


F658N (Ho-, [N II]) 


4 x 680 = 2720 





cluster populations in order to investigate the dependence of 
obtained LF parameters on maximum cluster mass, disruption 
time and cluster formation history. There, we also draw quali- 
tative conclusions about the cluster IMF and cluster disruption 
parameters of M5 1 . The variation of the bend luminosity with 
position on the disk is further investigated in §[5] and our con- 
clusions are summarized in § [6] 

2. Observations, photometry and sample selection 

2.1. Data 

We used the HST/ACS Hubble Heritage data of M51 (NGC 
5194) and its companion NGC 5195 in F435W, F555W, 
F814W (~B, V a nd I, respectively) an d F658N (Ha), see 
Table Q] (taken from lMutchler et alj J2005l) . proposal ID 10452, 
PI: S.V.W. Beckwith). The six pointings correspond to a total 
field of view of 430" x 610". This corresponds to 1 7.5 x 24.8 

covering 



kpc at a distance of 8.4 Mpc ( Fel dmeier et al 

f 



1993), 



a la rge fraction of the disk (more than 2 scale lengths of ~ 4 



kpc, Bec kman et all (119961) ) of M51 plus the region with NGC 
5195. This field is considerably larger than the field s of previ- 
ous studies dBastian et al.ll2005bL [oieles et al.ll2005l). 



to 



Mutchler et al 



Fo r details on reduction we refer 
( 2005b . The same data have also been used for our previ 
ous study of the lum inosity function of star clusters in M51 
(Gie es et al. 2006b) an d in a study on the radii of these clus- 
ters (IScheepmaker et al ■1 l2007h . 



For a detailed d escription of the data analysis we refer to 



Scheepmaker et alj (12007b . In view of the importance of the 
sample selection for this study a concise summary is given be- 
low. 



2.2. Source selection, photometry and radius 
measurements 

Selection of pointlike sources w as done with the SExtractor 
package ( Bertin & Arno uts 1996). Photometry was performed 
using the IRAF/DAOphot package. All magnitudes will be 
quoted in the Vega magnitude system. An aperture correction 
was applied, with the aid of artificial sources. The aperture cor- 
rections were all done for 3 pc sources. Since size and lumi- 
nosity ar e generally found to be only weakly correlated for star 
clusters Jzepf et aljll99St lLarsenll2004l iBastian etai1l2005bt 
Scheepmaker et al.ll2007l) . we do not expect the use of a con- 
stant aperture correction to introduce significant biases in our 
LF estimation. All clusters are affected by the same Galactic 



foregr ound extinction. We used the tables of Schle gel et al 



Table 2. 90 and 95% completeness limits for the three broad 
passbands. The same limits are used for different distance and 
background regions. They are all determined for a 3 pc cluster 
in a high background intensity region. 



Filter 




90% completeness 


95% completeness 


F435W (~ 


B) 


24.2 


23.6 


F555W (~ 


V) 


23.8 


23.0 


F814W (~ 


I) 


22.7 


22.2 



determinations one needs a wide range of broadband photome- 
try, in order to overcome the age/metallicity/extinction degen- 
eracy. Because we only have B, V,I photometry we are not able 
to clearly distinguish between the three effects. Therefore we 
can not correct for local extinction. 

We used the I SHAPE routine within the BAOlab package 
jLarsenll 19991 2004) to determine the effective radii (projected 
half light radii) of all point-like sources. Analytic cluster pro- 
files were convolved with an empirical PSF of the camera and 
fit to the sources. For high S/N (> 30) sources, sizes measured 
by ISHAPE on HST images are typically accurate to ~0.2 pix- 
els, so we adopt this a s the lower limit for an object we can 
recognize as resolved ( Larsenll2004 ). With ACS, at the distance 
of M51, this corresponds to 0.5 pc. We therefore take this as 
a lower limit. For a detaile d investigation of th e radii of this 
cluster population, see lScheepmaker et alJ (12007). 



(1998). For accurate local (i.e. in the M51 system) extinction 



2.3. Background regions 

Because the background intensity varies strongly over the im- 
age, we divided the image in three background levels, as indi- 
cated in Fig.Q] The image has been smoothed with a boxcar av- 
erage of 200 pixels. Two isophotes on this smoothed image are 
used as background limits (corresponding to 20.72 and 21.02 
mag arcsec~ 2 , respectively). We will use these three regions 
when investigating trends in the LF with background intensity. 
The high background region will also be referred to as being 
the spiral arms, whereas the low background will be called in- 
terarm region. 



2.4. Completeness 

Completeness tests were performed with the aid of artificial 
cluster experiments. In this study we make use of the complete- 
ness limits for 3 pc sources (see Table |2), determined in a high 
background region, to be sure that our sample is not affected 
by incompleteness due to high surface brightness. As a slope of 
a luminosity function might be sensitive to an incompleteness 
of ~ 10%, it is worth noting that using a higher completeness 
limit cuts off a big fract ion of the clusters. As can be found in 
Scheepmaker et alj d20071) . the completeness function becomes 
very flat above 90%. An estimate of the expected bias in de- 
termining the slope of the low luminosity end of the LF in the 
three filters can be found in Section[3] 
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Fig. 1. The contours that outline the three background intensity 
level regions, superimposed on the image in the F555W pass- 
band. The solid line encloses the highest background level and 
everything outside the dotted line is called low background. 
The area in between the two lines is referred to as being a tran- 
sition region, to have the other two regions clearly distinguish- 
able. 

2.5. Sample selection 

Finally, sources were selected from the complete sample (i.e. 
the sample of clusters, brighter than the 90% completeness 
limit) if: 

1 . the source is detected in F435W, F555W and F814W, 

2. the source was fit better with an extended cluster profile 
than with a delta function, meaning that the xl of the ex- 
tended profile is lower than that of a pure PSF fit, in the fil- 
ter of interest (this results in slightly different cluster sam- 
ples in different filters), 

3. the cluster has a radius (in the case of an elliptical cluster 
profile, the semi-major axis is converted to the radius of a 
circle with the same area) of at least 0.5 pc, 

4. the nearest neighboring point source is at least 5 pixels 
away to reduce the contamination. 

The resulting sample contains 7698, 6846 and 2539 clusters 
for F435W, F555W and F814W, respectively. To be sure we 



use a reliable sample of clusters we visually inspected the 
sample brighter than the determined completeness limits that 
pass the criteria above. For a description of the criteria to 
remove clusters from the sample and a graphical represen- 
tation of clusters t h at are removed by visual inspection, see 



Scheepmaker et al.l (120071) . Note that the sample used in this 



study still co ntains more sources than t he sample in the study 
described by IScheepmaker et al.1 (120071) . since there they need 
to be sure that they have a reliable estimate for the radius of the 
source, whereas in this study the only constraint is that we need 
to be sure that the source is extended. 

3. The luminosity function of clusters in M51 

Following the source selection as described in Section [2] 
we obtain the luminosity function of the star cluster 
candidates in M51. The distributions are fit with sin- 
gle and double power laws. In Appendix [A] we com- 
pare several methods for determining the shape of the LF. 
Following the results we will use a method that is de- 
sc ribed and tested bv iMafz Apellaniz & Ubeda (2005), based 
on lD'Agostino & Stephensl (ll986). In short, it uses bins that are 



variable in width, such that every bin contains an equal number 
of clusters. This avoids differences in statistical weights when 
fitting the distribution function, which would bias the obtained 
slope of the power law to bins with high statistical weight. The 
power law LF will be defined as 



dN 



dL A 

or, equivalently 



dN 



10 s ' 



A/, 



(1) 



(2) 



where L A (M A ) is the luminosity (magnitude), N is the number 
of clusters and a and /3 are related as a = 2.5 • /? + 1. 

Equation|2]is applied in fits for single power laws. For dou- 
ble power laws we apply 



dN 



for M A > M bend 
for M A < M bend 



(3) 



in which both power law slopes, the location of the bend 
(A^bend) and the normalization are free parameters. Note that 
only the normalization of one power law is a free parameter, 
the other one is then set by the location of the bend and that 
normalization. The total number of free parameters for the sin- 
gle power law is 2, while for the double power law fits there 
are 4 free parameters. The number of bins, and therefore the 
number of degrees of freedom, is different for all samples. 
The number of bins (M,ins) is related to the number of objects 

2/5 

in the sample (N r ) by Nh™ = 2 ■ N c +50 (adapted from 
D'Agostino & Stephens!! 19861) . 

The fit is performed on the whole range from the 90% com- 
pleteness limit to the brightest cluster in the sample. As we 
use the 90% completeness limits, we will here estimate by 
how much we might underestimate the slope of the low lu- 
minosity end due to incompleteness. By correcting the LF at 
the 90 and 95% limits for the incompleteness the LF would 
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become steeper. If we were to correct the value of the LF at 
the 90 and 95% completeness limits, the change in /3 would 
be A/? = 2.5/AM A ■ (log(l/0.9) - log( 1/0.95), in which AM A 
is the difference between the magnitudes of the 90 and 95% 
completeness limits per filter. As Aor = 2.5 ■ A/3, we obtain 
the following estimates for the underestimate of the low lumi- 
nosity end slope of the LF, for F435W, F555W and F814W 
respectively: 0.098, 0.073 and 0.12. Note that these estimates 
are upper limits, as the 95% completeness limits always falls 
well within the faint end of the LF, as will become clear in the 
following section. The underestimate only holds for the range 
in magnitudes between the 90 and 95% completeness limits, 
brightward of the 95% limit the underestimate is less severe. 



3.1. Luminosity functions of the whole sample 

For F435W, F555W and F814W the LFs are shown in the top 
panels of Fig.[2l[3]and[4]respectively. They are fit with a double 
power law distribution function. This fits the data better th an a 
single power law, as already found bv lGieles et alj d2006bl) . 

To confirm their conclusion, the fits are also performed with 
a single power law distribution function, using the fit method 
as described in the Appendix. In order to compare the fits, a 
comparison of the (reduced) xl is made. There are two more 
parameters: an extra slope and a bend magnitude. Note how- 
ever, that these parameters are not independent, as the slope of 
one of the power laws, together with the location of the bend, 
define at least one point in the other power law (they meet at 
the bend location). The difference in fit quality is indicated by 
both reduced^ 2 (jfy in the final two columns, together with the 
number of degrees of freedom. 

In the upper panels of Fig. [2j [3] and [4] we also show the LF 
for all clusters with an effective radius bigger than 2 pc using 
dot-dashed lines. The LF parameters agree within lcr with the 
LF for all clusters with r e ff > 0.5 pc. This reassures that the 
sample is not contaminated with bright stars. 

For an overview of the fit results, see Table[3] Faintward of 
the bend magnitude the slope is shallower than brightward of 
the bend as expected f r om the artificial cluster populations pre- 
sented in Gie les et alJ (12 006a) and i n agreement with pr evious 
observations ( Whit more et all 1 9991: iGieles et al.ll2006bh . 



Hwang & Led (120081) used the same data set, but different 
cluster identification criteria, to determine the LF of the same 
population of clusters. Comparing th eir single power law fi t and 
the double power law fit as found bv lGieles et al. (2006b), they 
conclude that a single power law gives a better fit, by compar- 
ing X 2 f° r t> om fits on the F555W image. iHwang & L ee (2008) 
found for F555W on the interval -10 < Mfsssw < -8, a single 
power law with power law index -2.59 + 0.03. This is consis- 
tent with our high luminosity end slope of -2.52 + 0.06, which 
is found for Mp-^w < -8.23. The bend luminosity found in 
Gieles et al] d2006bl) is 0.7 mag brighter in F555W than the 
bend found in this work. The reasons are a correction for aver- 
age extinction (0.3 mag in F555W), a different bin ning method 




5 E~ 



M,= -6.70± 0.14 
a,= 1.73 ± 0.07 



2 r 

:5.5 < D < 9 kpc 



E_ 




and slig htly different cluster selection criteria in Gie les et al 
J2006bb . 



Fig. 2. The LF (number of clusters divided by the bin width) 
of the identified clusters in M51 in Mfa^sw- The double power 
law fits are performed on all of these clusters. Only clusters 
with a galactocentric distance < 9 kpc are in this sample. The 
number of clusters in the sample is indicated in all panels. The 
top panel is the whole sample. The lower three are samples, in 
number more or less equally divided, at different galactocentric 
radii. Both slopes and the magnitude of the bend are indicated 
(vertical dashed line). The dot-dashed line in the upper panel is 
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Table 3. Fit results of the whole sample in all three passbands. Every cluster brighter than the completeness limit in the filter 
under consideration is taken into account in the fits. Single power law as well as double power laws are fit, to judge the statistical 
validity of the bend. The first column is the passband, the second the number of clusters within the fit range. Column three 
contains the slope of the single power law fit, whereas the fourth, fifth, sixth and seventh column contain both the slopes and the 
location of the bend of the double power law respectively. The final two columns show the reduced^ 2 for single (xl s ) and double 
(xlj) power law fits, respectively, with the number of degrees of freedom indicated between brackets. All errors are lcr. 







Single 


power law 


Double power 


law 












Passband 


N 


a 






a 2 




M bend 




xl (#d.o.f.) 


x\ 4 (#d.o.f.) 


F435W 


7698 


1.98 ± 


0.01 


1.70 ±0.02 


2.34 


±0.04 


-7.49 ± 


0.09 


3.16(119) 


0.96(117) 


F555W 


6846 


2.07 ± 


0.01 


1.86 ±0.02 


2.51 


±0.06 


-8.23 ± 


0.12 


2.41 (116) 


0.99(114) 


FS14W 


2539 


2.28 ± 


0.02 


1.92 ± 0.01 


2.61 


±0.06 


-8.32 ± 


0.09 


1.50 (94) 


0.88 (92) 



Table 4. Fit results of the subsets in galactocentric distance in all three passbands. For every distance bin, the best fit faint- and 
bright end slope (a\ and ai, respectively) as well as the magnitude of the bend are indicated. All errors are lcr. The reduced^ 2 
values for both single and double power law fits (and the number of degrees of freedom) are given in the last column. 



Passband 


D G 


N 






Q-2 




M bend 


xi, s (#d.o.f.) 


Xij (#d-o.f.) 


F435W 


D a < 3 kpc 


1696 


1.30 


±0.04 


2.43 ± 0.09 


-8.09 ±0.10 


3.34 (87) 


1.00 (85) 




3 < Do < 5.5 kpc 


2975 


1.65 


±0.00 


2.41 


±0.03 


-7.16 ± 0.05 


2.11 (97) 


0.87 (95) 




5.5 < Dq < 9 kpc 


2669 


1.73 


± 0.07 


2.29 


± 0.05 


-6.70 ±0.14 


1.39 (94) 


0.94 (92) 


F555W 


D a < 3 kpc 


1607 


1.38 


± 0.00 


2.48 


± 0.02 


-8.28 ± 0.02 


2.87 (86) 


0.96 (84) 




3 < D G < 5.5 kpc 


2608 


1.82 


±0.04 


2.46 


±0.07 


-7.40 ±0.14 


1.75 (94) 


1.09 (92) 




5.5 < Dq < 9 kpc 


2286 


2.03 


±0.01 


2.38 


±0.08 


-7.51 ± 0.22 


1.19(92) 


1.08 (90) 


F814W 


D a < 3 kpc 


847 


1.53 


± 0.09 


2.65 


± 0.14 


-8.56 ± 0.14 


1.57 (77) 


0.84 (75) 




3 < Dq < 5.5 kpc 


883 


1.88 


± 0.18 


2.66 


± 0.11 


-7.80 ±0.18 


1.12 (78) 


0.93 (76) 




5.5 < D a < 9 kpc 


711 


1.97 


± 0.04 


2.44 


±0.09 


-7.48 ± 0.22 


0.81 (75) 


0.82 (73) 



Table 5. Fit results of the subsets in the high and low background regions in all three passbands. For both backgrounds, the best 
fit faint- and bright end slope (a\ and Qf2, respectively) as well as the magnitude of the bend are indicated. The final two columns 
show the reduced^ 2 for single Of 2 ,) and double Okfjjj) power law fits, respectively, with the number of degrees of freedom 
indicated between brackets. All errors are lcr. 



Passband 


Background 


N 


ai 




a 2 




M bend 


Xi, s (#d-o.f.) 


Xi, (#d.o.f.) 


F435W 


Low 


2837 


1.99 


±0.06 


2.49 


±0.07 


-6.91 ± 0.17 


5.28 (96) 


1.02 (94) 




High 


2771 


1.10 


±0.04 


2.24 


±0.05 


-7.41 ± 0.07 


1.19(95) 


0.83 (93) 


F555W 


Low 


2378 


2.19 


±0.08 


2.54 


±0.08 


-7.43 ± 0.29 


4.62 (92) 


1.00 (90) 




High 


2606 


1.38 


±0.07 


2.41 


±0.07 


-8.00 ± 0.08 


1.21 (94) 


1.11 (92) 


F814W 


Low 


623 


1.98 


±0.33 


2.61 


±0.10 


-7.30 ± 0.22 


1.57 (74) 


0.86 (72) 




High 


1333 


1.60 


±0.09 


2.50 


±0.08 


-8.20 ±0.13 


0.75 (83) 


0.73 (81) 



3.2. Variations with galactocentric distance 

The total sample is large enough to cut it in three, more or less 
equally sized (in number of clusters), galactocentric distance 
'bins', without losing too much statistical confidence. The lim- 
its are chosen at 3 and 5.5 kpc. For the distance dependencies 
of LF parameters, only clusters within a circle with galacto- 
centric distance smaller than 9 kpc are used, in order to be sure 
that all the clusters belong to M51 and not to its companion 
NGC5195. In the resulting three regions the luminosity func- 
tions are shown in Fig. [21 Fig.[3]and Fig.|4]for F435W, F555W 
and F814W respectively. The results of the double power law 
fits are shown in Table [4] 

We find that the location of the bend shifts to fainter magni- 
tudes for greater galactocentric distance for all filters. The bend 
magnitudes in the different filters get brighter for redder filters. 



Furthermore, in all three filters, the slope of the faint end of 
the LF gets steeper for subsets further out in the galaxy. In the 
D < 3 kpc subsamples it is even clearer that a double power 
law fits the distribution better, as the difference between slopes 
is biggest here. 

3.3. Variations with background intensity 

The background intensity regions, as described in Sect. 12.31 
also divide the sample in sub-samples, in which we can look 
for variations in the shape of the LF. The results are shown in 
Fig. [5] for the F435W passband. Because the results are very 
similar for different filters, we do not show the LF for other fil- 
ters, but we do list all results in Table [5] Only the high and low 
background regions are compared. 

The location of the bend is about half a magnitude brighter 
in higher background regions, although uncertainties are big- 
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more strongly in subsamples for which the faint- and bright- 
ger than for the galactocentric distance dependence. The faint- end slopes are strongly different). 

end slope depends quite strongly on background intensity. It If clusters are formed in the spiral arms, a younger popula- 
is much shallower in high background regions, again showing tion of clusters would be expected in the high background in- 
clearly the difference in fit quality between single and dou- tensity region. This is, for this sample of clusters, confirmed by 



ble power law descriptions (which of course manifests itself IScheepmaker et al.l d2007|) . We would therefore expect to find 
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Fig. 5. Same as Fig. [2] but now with varying background inten- 
sity, comparing spiral arms with interarm regions. 



in the high background region a faint-end slope of the LF, more 
closely resembling the assumed cluster IMF slope of -2. Th is 



slope is found in earlier stu dies (IBik et al.l l2003; Ba stian et al 
2005b; Gieles et al. 2006a). This is not what we observe. A 
possible explanation is the blending of faint sources. If one 
source is very close to another, and those are together selected 
as one source, then the light is added. Thereby, two fainter clus- 
ter disappear from the sample and one brighter cluster appears, 
which flattens the LF. It is difficult to estimate the detailed ef- 
fect of blending on the LF, since the youngest c lusters tend to 
be found within complexes (IBastian et al . 2005a) and blending 
will be a strong function of cluster age, luminosity and environ- 
ment. Although we have carried out completeness experiments 
in regions of high background, these might still be optimistic 
for the most extreme cases of crowding. 

The trends with background intensity are correlated with 
the trends with galactocentric distance, because the high back- 
ground regions are on average closer to the center than the 



regions of low background. It is hard to disentangle the two 
trends and judge which is the most fundamental of the two. We 
conclude that we find a significant dependence of the faint-end 
slope on whether the population of clusters is in a high or low 
background region and that the bend in the LF occurs roughly 
half a magnitude brighter in the high background intensity re- 
gion. 

4. Shaping the LF of the M51 cluster population 

To get a handle on the variation of the LF parameters across 
the disk, we will discuss the variation of the LF parameters 
of artificial cluster models. These models compute the inte- 
grated properties of a cluster population, such as the LF, age 
distribution and mass distribution, based on an assumed clus- 
ter formation history, initial cluster mass function and cluster 
disruption scenario, which are all variable. Fitting the observed 
LF parameters in the disk of M51 in detail requires fitting a 
large parameter space in cluster IMF, cluster formation history 
and disruption parameters as a function of time and position in 
the disk, especially for interacting galaxies, where star forma- 
tion histories and local environmental conditions are generally 
complex. 

Here, we will investigate a range in cluster formation histo- 
ries, disruption times and maximum masses in a Monte Carlo 
simulation. 



4.1. Artificial cluster populations 

If a certain cluster formation history and cluster IMF are as- 
sumed, one obtains a full history of cluster formation, which 
provides one with an age and initial mass for every clus- 
ter in the system. Together with a description for the clus- 
ter disruption (mass as a function of age for a cluster with 
certain initial mass), this gives at all moments the present 
day mass and age of a cluster. For given metallicity and stel- 
lar initial mass function, also the luminosity in every photo- 
metric passband is available. These come from stellar popu - 



lation synthesis models. The models of iGieles et alJ d2005h , 
that are used here, use the GALEV models (Schu lz et al. 2002; 
Anders & Fritze-v. Alvenslebenll2003l) . 



4.1.1. Cluster formation histories 

We construct artificial populations for different cluster forma- 
tion histories. The underlying smooth cluster formation history 
is taken to be a power law in time: 



dN 
~d7 



(4) 



where the parameter p is taken zero (constant cluster forma- 
tion in time), positive (p = 0.5, for a cluster formation rate in- 
creasing with time, i.e. higher now than in the past) or negative 
(p = -0.5, for a declining cluster formation with time). The 
ages are drawn between a minimum and maximum age of 10 7 
and 10 1() yr, respectively. On top of all these smooth formation 
histories we add an extra burst of cluster formation, forming an 
extra factor of (no burst) or 0. 1 (forming 10% of the mass of 
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the smooth population in the burst on top of the smooth forma- 
tion). This value of 0.1 is chosen to be an extreme case, as M51 
and its companion are not really merging, but rather a system 
in a close encounter. The age of the clust ers formed in the burst 
is taken to be 3 ■ 10 7 yr, the age at which iBastian et al] (]2005b) 
found a peak in the cluster formation history of this galaxy. 

We assume all cluster populations to form with a cluster ini- 
tial mass function that is described by a power law with index 
-2, truncated at the maximum mass (see below). 

4.1 .2. Maximum masses and disruption times 

Disruption times are taken log(f4/yr) = [oo, 9, 8] to investigate 
the difference between no disruption, moderate disruption and 
strong disruption. t$ is the total disruption time of a cluster with 
an initial mass of M,- = 1O 4 M . For other masses the disruption 
time is given by 



fdisCMi) = ? 4 (M,/lO 4 M y 



(5) 



with 



0.62 



(Boutl oukos & Lamersl I2003 
Baumgardt & Makinol I2003I) . For the mass of a cluster as 



a function of age we ta ke the gradual dissolution prescription 
of lLamers et alJ (120051) . without taking into account the mass 
lost due to stellar evolution (because the population synthesis 
models use initial stellar mass as input). The maximum masses 
are varied between 10 425 M and 10 6 M Q taking steps of 
log(M) = 0.25. The lower limit is chosen such that the clusters 
are still within the observable range, the upper limit goes up to 
the point where the mass function is not any longer sampled 
up to the maximum mass. This last point is disruption time 
dependent as will become clear in the following sections. 

4.1 .3. Analysis of the model population 

For all cluster masses a maximum lifetime is determined 
from their initial mass and their corresponding disruption time 
(Eq. O. For all surviving clusters we derive the magnitude, by 
scaling th e magnitude fro m the GALEV stellar population mod- 
els with a lKroupa (2001) stellar IMF for solar metallicity to the 
current mass. Because we do not have any information about 
internal extinction, we do not apply any attenuation on the clus- 
ter magnitudes. This may shift the bend magnitude of the LF 
to brighter magnitudes, by a few tenths of magnitudes, com- 
pared to the observations. Differential extinction may change 
the slopes of the power laws as well. 

Because the resulting LF parameters from the fits will de- 
pend on the number of clusters observed (which determines 
how far the LF is sampled at the high luminosity end) we follow 
a very big population of clusters. We make sure that the popu- 
lation consists of more clusters than we have observed. When 
fitting the LF in the same way as the data, we only take clus- 
ters brighter than our 90% completeness limits into account. In 
order to have the same statistics we now take a random sample 
of clusters from this population, exactly equal to the number of 
clusters observed in the specific filter under consideration. We 
repeat every run 20 times with a different random number seed. 

We apply the exact same analyses as for the observed data, 
fitting both single and double power laws to the distribution 
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Fig. 6. The bend magnitude as a function of maximum mass, 
for disruption times log(?4/yr) = [oo, 9, 8]. The cluster forma- 
tion rate is taken constant in time and there is no extra burst 
included. The three different line styles correspond to differ- 
ent disruption times as stated in the legend. The thick line is a 
Monte Carlo realization with the number of clusters in sample 
equal to the observed number of clusters in the same filter. The 
thin lines (without symbols) are for the same disruption time, 
now with the number of clusters in the sample scaled with the 
maximum mass. The thin long-dashed line is the line for which 
A^bend 00 -2.5 log(M max ), the relation expected on simple theo- 
retical grounds, as explained in the text. 



of cluster luminosities. We keep track of the difference in fit 
quality and will not include populations that are fit better with 
a single power law (determined from xl) or that have a better 
fit with a double power law, but with the bend at the edges 
of the range of luminosities in the sample (essentially also a 
single power law). We take the median bend magnitude of all 
the populations that satisfy these criteria. Furthermore, we only 
include realizations which sample the mass function up to the 
upper mass limit. If not sampled that far, one does not expect a 
bend in the LF and the bend that might be found does not mean 
much. Note that just a few runs are thrown out (5 out of 20 at 
most, in 90% of the t\ - M max combinations). 



4.2. Resulting LF parameters 

We will only show results for the F435W filter, as here the 
statistics of the observations are best due to the longer integra- 
tion time. Results are not qualitatively different in other filters, 
although the exact value of the slopes of the LF and the bend 
magnitude are different. 

In Figure [6] the resulting bend magnitude as a function of 
maximum cluster mass is shown for the three different disrup- 
tion times, for a constant cluster formation rate, without any 
extra burst. It can be clearly seen that, in general, the bend 
shifts to brighter magnitudes for higher maximum masses. At 
high masses, the relation is cut-off, because the mass function 
is no longer sampled up to the maximum mass, when using 
the same number of clusters as observed. When repeating the 
experiment where we scale the number of clusters in the sam- 
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pie linearly with the maximum mass (thin lines in the figure), 
the upward trend is continued. In this case the mass function is 
always sampled to well above the upper mass limit. 



4.2.1 . Explaining the dependence of the bend 
magnitude on maximum mass 

For a cluster sample with a uniform age distribution and no dis- 
ruption, the bend luminosity is that of a cluster of the maximum 
allowed mass, with an age equal to the age of the cluster sample 
dGieles et al.ll2.005b . If disruption effects and non-uniform age 
distributions are allowed the relation between bend luminosity 
and maximum mass becomes more complicated, as discussed 
below. 

For no disruption, we can see in Figure [6] (dotted line) 
that, if the maximum mass is considerably higher than the 
lower mass limit of visible clusters (depending on age be- 
tween 10 3 and 10~ 45 ), the relation between maximum mass 
and bend magnitude is given by Mb enc i K -2.5 log(M max ), 
as expected from simple arguments: luminosity depends lin- 
early on mass for fixed age and metallicity, and the magnitude 
M,\ oc -2.5 \og{L,\). This relation is overplotted on Figs. [6] and 
[7]with the long dashed lines. 



4.2.2. The effect of cluster disruption on LF 
parameters 

The main influence of cluster disruption is seen in the slope 
of the low-luminosity end of the LF. This is because the clus- 
ter disruption time depends on its mass, with low mass clusters 
disrupting faster than more massive clusters. Disrupting a clus- 
ter means moving it to a lower mass bin, so mass dependent 
disruption means that clusters are moved to a bin belonging 
to lower mass at an ever decreasing rate for increasing mass. 
Cluster disruption therefore flattens the LF and stronger dis- 
ruption will lead to stronger flattening. 

The bend location is also affected by the disruption time. 
The maximum mass occurring in a log(age) bin, now is not the 
maximum initial mass, but the maximum initial mass evolved 
in time due to disruption. In case of an untruncated mass func- 
tion, this would just make the relation between the number of 
objects in a log(age) bin grow slower than linear with age, re- 
sulting in a shallower LF slope, than the slope of the cluster 
initial mass function. When adding a maximum mass, the high- 
est possible mass is also a function of age. Stronger disruption 
lowers the mass of a cluster more during the lifetime of the 
galaxy. Therefore, for the same maximum mass the bend be- 
comes fainter at the high maximum mass end of the Figure 
©. For low maximum masses, a cluster of the maximum mass 
can be completely disrupted in the lifetime of a galaxy. In that 
case, the bend would occur at the luminosity of the most mas- 
sive surviving cluster at a relatively young age. In that case, the 
bend is brighter than it would be without disruption, see Fig. [6] 
Summarizing this section we conclude that disruption makes 
the relation between the maximum mass and bend magnitude 
flatter than it would be without disruption. 
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Fig. 7. Same as Figure [6] but now for three different cluster 
formation histories, upper panel: an increasing cluster forma- 
tion rate with time (dN/dt oc f - 5 ) without burst, middle panel: 
a decreasing cluster formation rate (dN/dt oc f -0 - 5 ), without 
burst and lower panel: a constant cluster formation rate, as in 
Figure |6]but with an extra burst, forming 10% of the stars with 
an age of 3 • 10 7 yr. 



4.2.3. The influence of the cluster formation history 

The factor that is probably worst constrained by the data is 
the cluster formation history. Here we will show that a differ- 
ence in cluster formation, even in simple cases, already makes 
the obtained bend magnitude as a function of maximum mass 
very degenerate with disruption and maximum mass. We sim- 
ulate a population of clusters with the power law cluster for- 
mation rate as a function of time given by Equation |4] with 
p = 0.5 and p = -0.5, to investigate the resulting bend magni- 
tudes for cluster formation rates increasing and decreasing with 
time, respectively. We also investigate the effect of a starburst 
on top of a constant cluster formation history, as described in 
Section 14.1.11 The resulting maximum mass-bend magnitude 
relations are shown in Figure [7] 

The behavior of the bend magnitude with maximum mass 
for a decreasing cluster formation rate with time (i.e. lower now 
than in the past) is very similar to that of a constant star forma- 
tion history (Fig. [6}, as is shown in the middle panel of Fig. [7] 
As the bend occurs at the magnitude of the highest mass clus- 
ter at the age of the population, this is expected, since in both 
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cases clusters with relatively old ages receive strong weight in 
the combined LF of the whole sample. Therefore the bend mag- 
nitude, which is the magnitude below which clusters of all ages 
contribute to the LF, remains unchanged. 

For a cluster formation rate which is increasing with time, 
the oldest clusters receive relatively little weight in the com- 
bined LF, and therefore the bend does not show up clearly at the 
expected magnitude. Instead, the LF is dominated by the com- 
bination of cluster disruption and the cluster formation history. 
For stronger disruption, the LF becomes more weighted to- 
wards younger clusters, so for stronger disruption the obtained 
bend is brighter. 

An extra burst, as depicted in the lower panel of Figure [7] 
results in very irregular relations between bend magnitude and 
maximum mass. It depends strongly on the fraction of stars 
formed in the burst and the age of the burst. Therefore, with- 
out a very detailed knowledge about the burst, constraining the 
mass function or disruption is nearly impossible. 

4.2.4. The wavelength dependence of the bend 

The bend luminosity is expected to be higher in redder filters, 
as the shape of a clusters spectral energy distribution is dom- 
inated by red stars for clusters with the age of the population 
(i.e. old clusters). In red filters there is least evolution in the lu- 
minosity over the lifetime of a cluster. This results in a steeper 
bright-end slope for the red filter than for the blue filter. 

As described in iGielesI (12008b . the bright end slope of the 
LF can be predicted and depends only on the fading of a clus- 
ter in a specific filter. If the luminosity of a cluster in a fil- 
ter as a function of time can be approximated by a power 
law, as L(A, t) oc r^, then the bright end LF slo pe should be 
dN/dL^ un oc Lj 1 Taking values for £ from iGieles et al 



(2007b) for clusters younger than 3 Gyr), we obtain the 
following predicted high luminosity end slopes for F435W, 
F555W and F814W: -2.28, -2.52 and -2.85, respectively. 

4.3. The cluster population of M51 

On the basis of a comparison between, on the one hand, the 
qualitative explanation of LF parameters as a function of clus- 
ter IMF and cluster disruption given above, and on the other 
hand the observed trends listed in Tables [3] E]and|5]the follow- 
ing qualitative interpretation of the observed LF parameters can 
be given: 

1 . The disruption time for star clusters varies with galacto- 
centric distance, such that disruption is strongest in the 
central regions of the galaxy. This conclusion can be drawn 
from the shallower faint-end side of the LF in the central re- 
gions of the galaxy. Concluding on the dependence of dis- 
ruption time on spiral structure is hampered by the contra- 
diction between expected and observed dependence of the 
faint-end slope of the LF on spiral structure. Besides, there 
is a degeneracy between variations as a function of galacto- 
centric distance and background intensity that arises from 
the fact that high background intensity regions preferen- 
tially lie closer to the center of the galaxy than high back- 



ground intensity regions. Selection effects, though, may be 
important. 

2. Translating trends in bend magnitude with radius and back- 
ground level to a trend in terms of maximum cluster mass is 
non-trivial and requires detailed knowledge about the for- 
mation and disruption histories of the cluster populations 
involved. We tentatively suggest that the observed trend 
of decreasing bend magnitude with galactocentric distance 
(and lower background level) may be indicative of a corre- 
sponding decrease in maximum cluster mass. Assuming a 
constant upper mass limit and simple cluster formation his- 
tory would require large variations in the disruption time 
across the disk. 

The conclusion on differences in disruption time between 
spiral arms and interarm regions might not hold if the age 
distribution in spiral arms is different from the age distri- 
bution in interarm region s. This is found to be the case by 
Scheepmaker et alj J2.007I) . The variation of the bend luminos- 
ity between arm and interarm regions might be explained by 
extinction, which is probably stronger in the spiral arm than in 
the interarm regions. With this limited set of broad band filters 
we can not test this hypothesis. 



5. The variation of the bend luminosity across the 
disk 

Here, we investigate the variation of the luminosity of the bend 
over the disk. At fixed age, the luminosity of a cluster scales lin- 
early with its mass, and therefore the variation of the luminos- 
ity of the bend should follow the variation in maximum mass 
across the disk. 

As luminosity scales linearly with mass and magnitudes 
logarithmically, we will fit the quantity L = 10~°- 4 ' Mw as a 
function of galactocentric distance. Fitting an exponential to 
the three points in galactocentric distance (the mean galacto- 
centric distance of all clusters in the sample; note that this 
is different for the three filters, as the samples are not the 
same), versus log Lbend (the logarithm of luminosity in units of 
the filter dependent zero point luminosity), gives scale lengths 
of 3.4+0.4, 4.5+0.7 and 4.7+1.0 kpc for F435W, F555W and 
F814W, respectively. See Fig. [8] In these fits, the errors on 
the bend magnitudes are taken into account. The error on 
mean galactocentric radii for all sub-samples are assumed to 
be negligible (the error can be estimated to be a few pixels 
for every cluster, which results in relative errors much smaller 
than the relative errors on the bend magnitude). When fitting 
power laws to the points, the indices (p) of the power laws 
(log(L bend ) cc p ■ \og{D ga ,)) for the fits through F435W, F555W 
and F814W are -1.0+0.1, -0.8+0.1 and -0.7+0.1 respectively. 
These fits are also shown in Fig. [8] 

Because we have only three points per filter, and because 
of the lack of detailed modelling of the cluster population (and 
therefore the lack of precise determinations of the maximum 
mass), we do not want to discriminate between a power law 
dependence and an exponential decline of the bend luminosity 
with galactocentric distance. 
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Fig. 8. For all three filters, -OA-Mbend = log(Lt, en d/ioX with Lq 
the filter dependent zero-point luminosity is plotted against the 
galactocentric distance. Exponential fits are shown for all three 
filters, as well as power law fits. The error bars are lcr errors. 

6. Conclusions 

Using high resolution, deep HST/ACS imaging, we have con- 
firmed the fact that the luminosity function of star clusters in 
M51 is better described by a double power law distribution 
function than by a single power law. Although any sub-sample 
of clusters is fit only slightly better by a double power law than 
by a single power law, the fact that it holds for all sub-samples 
strengthens the belief that the LF is indeed fitted better with 
a double power law than with a single power law. The fits on 
small galactocentric distance bins and high background regions 
clearly show that a single power law is insufficient in describing 
the LF. Also, the fact that for every subsample the luminosity of 
the break is higher in filters with longer effective wavelength, 
as expected due to evolutionary fading of the clusters strength- 
ens the conclusion that the bend in the LF is connected to a 
maximum mass (by comparison with the models). The bright 
end slope of the LF is in all filters c onsistent with the slope 
expected from analytical arguments of lGielesI (|2008). 

The interpretation of the double power law distribution 
function is guided by the artificial population models and leads 
to the following conclusions: 

1 . The bend in the LF occurs at brighter magnitudes, closer to 
the center of the galaxy. This may imply a higher possible 
maximum mass closer to the center of the galaxy. This con- 
clusion is guided by the Monte Carlo simulations, which 
show that for given cluster formation history and disrup- 
tion time, the bend magnitude is brighter for higher maxi- 
mum masses. Varying the cluster formation history, and/or 
cluster disruption across the disk makes this inference less 
straightforward. However, assuming a constant maximum 
mass would require large variations in the disruption time 
as a function of galactocentric distance. 

2. The magnitude of the bend in the LF varies slowly with 
background intensity. This variation may be solely due to 
the difference in average galactocentric distance between 
both background intensity regions. 



3. The other deduction from the LF parameters concerns dis- 
ruption parameters, which also depend on galactocentric 
distance. A higher surrounding density will be more de- 
structive for a cluster, as encounters with massive objects, 
clouds and clusters, will be more frequent. Also, the tidal 
field of the galaxy is decreasing outwards. Due to the de- 
generacy of a decreasing upper mass limit and stronger dis- 
ruption to flatten the faint end of the LF, the dependence 
of the disruption time on galactocentric distance can not be 
constrained. This degeneracy can be lifted if the ages and 
masses of the clusters can be measured. 

4. The slope of the faint end side of the LF is shallower, closer 
to the center of the galaxy and supports the suggestion that 
star clusters in the center are more easily disrupted. The 
higher density, stronger tidal field and larger relative ve- 
locities between several massive components (clusters and 
clouds) make the di sruption of star clusters ea sier in those 
parts of the galaxy dGieles et alJl2006cll2007ah . 

5. The slope of the faint end side of the LF is shallower in 
high background regions. This is counter-intuitive, as one 
expects the younger population in the spiral arms to have a 
faint-end LF slope more closely resembling the cluster IMF 
slope compared to the population in the interarm regions. 
Selection effects may play an important role here. Clusters 
have a high number density, causing confusion, blending 
and the rejection of clusters from the sample as a result of 
a very nearby neighbor. 

Whether galactocentric distance or background intensity is 
the driving parameter for the variations of the faint-end slopes 
and bend magnitudes is unclear, due to the correlation of back- 
ground intensity with galactocentric distance. 
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Appendix A: Comparison of fit methods 

In order to fit a distribution function to the luminosity dis- 
tribution of the clusters we exploit a meth od described and 
tested by Mafz Apellaniz & Ubeda (120051) . who adopted it 
from D'Agostino & Stephens! d 19*86 ). In conventional fit meth- 
ods one uses equally sized bins (in either linear or logarith- 
mic scale) for the quantity in question, luminosity in our case. 
Because the distribution is described by a power law with neg- 
ative index, this always means that the bins at lower luminos- 
ity will contain more clusters and will therefore get a higher 
statistical weight in the fit (due to the Poissonian nature of 
the statistics in counting experiments). This results in an, on 
average, shallower slope from the fit routine than is actually 
pr esent in the data, as shown by t he Monte Carlo simulations 
of iMafz Apellaniz & Ubedal ( 20051) . In order to overcome this 
problem these authors propose to use bins of variable binsize, 
chosen such that every bin contains equally many clusters and 
therefore have the same absolute Poissonia n error. Besides the 
method of IMafz Apellaniz & Ubeda d2005l) we tested a cumu- 
lative distribution method, used by iRosolowskyl ( 2005 ) and a 
maximum likelihood method. 

In order to make a comparison between the four methods 
to fit a distribution function, w e follow a method similar to 
Mafz Apellaniz & Ubedal ( l2005h . For every method a sample of 
10.000 numbers is randomly drawn from a power law distribu- 
tion function with input slope -2.35, by linking a random num- 
ber from a uniform distribution to the cumulative probability 
density function, belonging to the power law distribution func- 
tion. This sample is fit by the four methods (described below) 
and the resulting slope is recorded (we fitted on the probability 
density function, such that there is only 1 free parameter: the 
slope). This process is repeated a 1000 times and the difference 
between the input slope (//input) and the slope that results from 
the fitting procedure (rjf\ t ) is plotted in Fig. IA.ll The methods: 



1. The 'conventional method' uses equally sized bins (either 
linearly or logarithmically, linearly is shown in the upper 
left panel of Fig. IA. lb . The problem of the difference in 
statistical weights makes the resulting histogram a function 
of the slope used for the distribution function (the steeper 
the function, the larger the difference in number of clusters 
in the low mass/luminosity bins and the greater the spread 
in Poissonian errors). In the case of a slope 77 = -2.35 this 
results in a large systematic offset as shown in Fig. IA.lK up- 
per left panel). 

2. Mafz Apellaniz & Ubedal (2005) use a method with vari- 
able bin sizes, such that every bin contains equally many 
objects. The fit is consequently performed on the width of 
the bin, instead of the height. Errors in the widths are very 
small if photometric errors are neglected (ordered masses 
or luminosities are very close together for large samples), 
such that the Poissonian error on the number of clusters 
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Fig.A.l. A comparison of four different methods to fit a distribution function to a data set. All frames are on the same scale, 
but the horizontal axis of the upper left panel is different. Histograms of the difference between the input slope and the fitted 
slope (?7fj t - 77 m p U t, with 77 m p U t = -2.35) are shown for four different methods to fit a power law distribution function. Upper 
left: equally sized, linear bins. Upper right: bins with an equal number of objects per bin. Lower left: cumulative mass function. 
Lower right: Maximum Likelihood. In all panels the vertical dashed line indicates the center of the Gaussian that is fit to the 
distribution of fit results. See text for details on the methods. 



in the bin (the square root of that same number) is always 
dominant and the same for every bin. Using this kind of fit- 
ting results in residual slopes as shown in Fig. IA. II (upper 
right panel). 



A method used, among others, by IRosolowskvl d2005l) . in- 
volves no binning and makes use of the equation 



N(m > m 



*Jm' 



N(m) dm 



(A.l) 



(the cumulative distribution). The left hand side is the num- 
ber of objects with a mass greater than or equal to m' 
and C is the normalization (from total number of clusters). 
While the fact of not having to deal with bins (and therefore 
Poissonian errors) is ideal, the computational time neces- 
sary to deal with this method becomes long (quadratically 
with the number of objects), and therefore the lower left 
panel of Fig. lA.ll is composed of samples with only 1000 
objects, ten times less than the other ones. In this study, the 



number of clusters is large and therefore we will not use 
this method. 

A reliable method that does not make use of bins is 



the m aximum likelihood method dBevington & Robinson 



2003). The likelihood that a certain guessed theoretical dis- 
tribution function underlies the data is maximized. Here, 
the 'guessed' distributions are power laws in which only 
the slope is changed in very small steps (0.005) from -0.35 
to -4.35. The likelihood function £, is given by 



(A.2) 



In here, / is the function describing the guessed distribu- 
tion, and Xj are the values of the quantity one wants to fit 
(e.g. the luminosity), for every object i in the sample. In the 
case of a power law distribution function for luminosities 
the likelihood function would be 



(A.3) 
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So for every input value of a the likelihood function will 
have a different value. The shape as a function of alpha will 
generally be gaussian. The peak of this gaussian defines the 
best fit for the power law slope, and the FWHM of the gaus- 
sian is the standard deviation. As the likelihood function 
generally consist of small numbers, it is common practice 
to maximize the logarithm of the likelihood function: 

log(£) = -or • J] log(L ; ) (A.4) 

i 

The maximum likelihood is found at offsets from the input 
value as shown in Fig. I A. 1 1 (lower right panel). Whereas 
this method in principle is the most reliable of the four, 
fitting several parameters simultaneously gives likelihood 
functions with several local minima, that may be compara- 
ble in significance. Whereas this would imply that a single 
good fit cannot be obtained, it is hard to automize the se- 
lection of the best fit. Therefore, and because the difference 
in fit quality with the method with variable bins is almost 
negligable, we choose not to use this method. 

For consistency we always use the same method for fitting the 
distribution function. We choose for the second method, be- 
cause we fit more than one parameters at once, in samples of 
thousands of clusters. 



